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ABSTRACT 

In  order  for  the  NIF  (National  Ignition  Facility)  to  achieve  its  most  challenging 
goal  of  fusion,  it  is  paramount  that  we  explore  new  laser  matter  interaction  regimes  and 
control  mechanisms.  This  thesis  explores  a  novel  interaction  regime  namely,  the 
enhancement  of  laser  absorption  when  an  overdense  plasma  is  irradiated  with  ultra- 
intense  laser  light.  Self-generated  magnetic  fields  of  up  to  1 09  gauss  have  been  predicted 
in  overdense  plasmas  irradiated  with  ultra-intense  laser  light  pulses.  We  find  that  laser 
light  absorption  can  be  significantly  enhanced  by  the  oscillation  of  electrons  across  these 
fields.  There  is  then  a  very  large  Lorentz  force,  which  can  strongly  accelerate  electrons 
into  the  plasma  and  can  lead  to  generation  of  harmonics  in  the  reflected  light.  We  also 
show  that  large  magnetic  fields  can  be  significantly  amplified  to  even  larger  values  by  the 
pressure  of  the  light  pulse.  Potential  applications  of  this  interaction  regime  include  the 
Fast  Igniter  approach  to  ICF  (Inertial  Confinement  Fusion),  the  creation  of  anti-matter, 
and  the  development  of  laser  radiography  to  name  a  few.  Zohar,  a  2-V*  D  PIC  (Particle  In 
Cell)  code  was  used  to  run  our  simulations.  Zohar  assumes  a  collisionless  plasma,  solves 
Maxwell's  equations  in  3D,  and  accounts  for  the  relativistic  effects  which  can  occur  in 
the  high-energy  regime  that  we  consider.  The  high  number  of  particles  simulated 
required  significant  computing  power,  which  was  obtained  by  remote  access  to  the  LC 
(Livermore  Computing)  East  computer,  which  provided  8  processors  and  8  GB  of  shared 
memory  at  440  MHz.  "Exceed"  software  was  required  to  graphically  display  our  results 
on  a  PC.  This  thesis  has  been  a  collaborative  effort  between  NPS  (Naval  Postgraduate 
School)  and  LLNL  (Lawrence  Livermore  National  Laboratories),  to  the  benefit  of  DoD. 
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I.  INTRODUCTION 

The  study  of  the  interaction  of  intense  laser  light  with  plasmas  serves  as  an 
excellent  introduction  to  the  field  of  plasma  physics  [Ref.  1],  Both  the  linear  and 
nonlinear  theories  of  plasma  waves,  instabilities  and  wave-particle  interactions  are 
important  for  understanding  laser  plasma  coupling.  Numerous  plasma  effects  have  been 
observed  in  laser  plasma  experiments,  and  many  challenging  problems  remain  to  be 
understood. 

A.  BASIC  PLASMA  CONCEPTS 

Plasma  is  a  medium  containing  many  charged  particles  governed  by 
electromagnetic  forces.  Plasma  is  a  quasineutral  gas  of  charged  particles,  which  possess 
collective  behavior.  Plasma  designates  the  fourth  state  of  matter,  which  exists  in  the  sun, 
stars,  and  space,  etc;  moreover,  most  of  the  Universe  is  made  up  of  plasma.  The  plasma 
medium  is  described  macroscopically  (on  a  large  scale)  by  its  temperature  and  density, 
and  changes  in  the  plasma  are  calculated  by  using  conservation  equations  such  as 
conservation  of  energy,  momentum  and  mass.  On  a  microscopic  (small)  scale,  a 
statistical  description  using  probabilities  is  required  to  describe  a  plasma  in  the 
calculation  of  the  positions  and  velocities  of  all  particles.  Due  to  their  mutual  collisions, 
the  plasma  charged  particles  emit  radiation  (electromagnetic  waves).  Moreover,  many 
different  waves  can  be  created  in  a  plasma  medium.  The  stabilities  and  instabilities  of 
these  waves  play  a  crucial  role  in  plasma  systems.  Plasma  plays  an  important  role  in 
explaining  our  universe  and,  as  we  shall  see,  scientists  believe  that  the  understanding  of 
plasma  is  the  key  issue  for  solving  the  energy  problems  of  mankind. 


A  plasma  is  basically  a  system  of  N  charges  that  are  coupled  to  one  another  via 
their  self-consistent  electric  and  magnetic  fields.  Following  the  evolution  of  these  N 
charges  in  principal  would  require  us  to  solve  6*N  coupled  equations  which  would  be 
tedious  at  best.  Fortunately,  great  simplification  is  possible  if  we  focus  our  attention  on 
the  collisionless  plasma  behavior. 

Let  us  assume  a  plasma  space  in  which  we  decompose  the  electric  field  into  two 
fields  (Ei  and  E2)  that  have  distinct  spatial  scales.  The  field  Ei  has  spatial  variations  on  a 
scale  length  much  less  than  the  electron  Debye  length  ADe .  The  field  of  an  individual 
charge  is  shielded  out  by  the  response  of  the  surrounding  charges  over  this  length.  Ei 
represents  the  rapidly  fluctuating  microfield  due  to  multiple  and  random  encounters 
(collisions)  among  the  discrete  charges.  In  contrast,  the  field  E2  varies  on  a  space  scale 
greater  than  or  comparable  to  XDe.  E2  represents  the  forces  due  to  the  average  or 
"collective"  motion  of  the  charges. 

We  thus  have  a  natural  separation  into  collisional  and  collective  behavior.  The 
collisional  behavior  becomes  negligible  when  the  number  of  electrons  in  the  sphere  with 
a  radius  equal  to  the  electron  Debye  length  becomes  very  large. 

1.    Concept  of  Temperature 

In  order  to  better  understand  the  complex  nature  of  plasma  behavior,  it  is 
necessary  to  expand  our  most  basic  notion  of  temperature  [Ref.  2].  We  note  that  a  gas  in 
thermal  equilibrium  has  particles  of  many  different  velocities  and  the  most  probable 
distribution  of  these  velocities  is  known  as  the  Maxwellian  distribution.  The  average 
kinetic  energy  (Eave)  in  three  dimensions  is: 

Eave  =  3/2KT  (1) 


where  K  is  Boltzmann's  constant,  K  =  1 .38  x  10"23  J/°K  and  T  is  the  temperature  in 
degrees  Kelvin.  The  general  result  here  is  that  Eave  is  J4  KT  per  degree  of  freedom. 
It  is  customary  in  plasma  physics  to  give  temperatures  in  units  of  energy.  To  avoid 
confusion  on  the  number  of  dimensions  involved,  it  is  not  Eave  but  the  energy 
corresponding  to  KT  that  is  used  to  denote  the  temperature.  For  KT  =  1  eV  =  1.6xl0"19  J, 
we  have  T  =  1 .6  x  1 0"19/l  .38x1 0"23  =  1 1 ,600,  thus  the  conversion  factor  is 
leV=l  1,600°  K. 

It  is  interesting  that  plasma  can  have  several  temperatures  at  the  same  time.  It 
often  happens  that  the  ions  and  electrons  have  separate  Maxwellian  distributions  with 
different  Tj  and  Te  (ion  and  electron  temperatures  respectively).  This  comes  about 
because  the  energy  exchange  between  colliding  electrons  is  larger  than  that  between  an 
ion  and  electron.    Since  colliding  electrons  are  the  same  size  and  mass,  appreciable 
kinetic  energy  exchanges  can  result  from  an  electron-electron  collision.  Now,  when  an 
electron-ion  collision  occurs,  we  recall  that  the  electron  is  only  1/3  700th  the  mass  of  an 
ion  and  therefore  a  finite  amount  of  time  (on  the  order  of  hundreds  of  microseconds  let's 
say)  is  required  for  a  couple  thousand  electrons  to  collide  with  an  individual  ion  to  result 
in  an  appreciable  kinetic  energy  exchange.    Bearing  this  in  mind,  we  can  see  how  each 
species  can  be  in  its  own  thermal  equilibrium  but  the  plasma  may  not  last  long  enough  for 
the  two  temperatures  to  equalize  (since  many  of  the  interesting  effects  of  the  experiments 
are  observed  in  nanoseconds). 

When  there  is  a  magnetic  field  B,  even  a  single  species,  say  ions,  can  have  two 
temperatures.  This  is  because  the  forces  acting  on  ions  along  B  are  different  from  those 
acting  perpendicular  to  B  (due  to  the  Lorentz  force).  The  components  of  velocity 


perpendicular  to  B  and  parallel  to  B  may  then  belong  to  different  Maxwellian 
distributions  with  temperatures  Tparaiiei  and  Tperpendicuiar. 

Interestingly  enough,  a  high  temperature  does  not  necessarily  mean  a  lot  of  heat 
[Ref.  2].  For  example,  the  electron  temperature  inside  a  fluorescent  light  bulb  is  about 
20,000°  K,  but  it  is  cool  to  the  touch.  This  is  because  heat  capacity  has  to  be  considered. 
The  heat  capacity  of  the  electrons  in  the  fluorescent  light  bulb  is  minimal  because  the 
electron  density  is  much  less  than  that  of  a  gas  at  atmospheric  pressure  and  the  amount  of 
heat  transferred  to  the  wall  by  electrons  striking  it  at  their  thermal  velocities  is  minimal. 
2.    Plasma  Instabilities 

Plasmas  behave  in  a  collective  way.  Different  types  of  waves  are  created  in  the 
plasma  due  to  the  collective  behavior  of  the  electrons  and  the  ions  [Ref.  2].  The  types  of 
waves  that  can  be  created  in  plasma  are  numerous  and  we  shall  mention  just  a  few  of  the 
most  fundamental  ones  to  the  study  of  plasma  behavior. 

Conditions  of  collective  motion  can  be  stable  or  unstable.  For  example,  suppose 
a  small  displacement  takes  place  in  plasma  confined  by  a  magnetic  field.  If  the 
collective  motion  behaves  in  a  stable  way,  the  displacement  will  be  restored  to  the 
original  motion.  However,  if  the  collective  motion  is  unstable,  any  small  displacement 
increases  rapidly  resulting  in  the  collective  motion  phenomenon  breaking  up. 

Numerous  instabilities  can  be  excited  by  an  intense  electromagnetic  wave  in  a 
plasma.  Many  possibilities  exist  even  without  dc  magnetic  fields  [Ref.  1].  These 
various  instabilities  can  be  thought  of  as  the  resonant  coupling  into  two  other  waves. 
Various  possibilities  and  their  names  are  shown  below. 

■     Ion  acoustic  decay  instability  o)0  -» to    +eot  at  n  =  iicr 


■  Raman  instability  co0  -»  cosc  +  cope  at  n  <  Va  iicr 

■  Brillouin  instability  cd0  ->  cosc  +cd1  at  n  <  ncr 

■  2  6^  instability  <y0  ->  6>pe  +  6>pe  at  n  =  lA  ncr 

Here,  co0  is  the  incident  light  frequency,  cosc  is  the  frequency  of  a  scattered  light  wave 

cope  is  the  frequency  of  a  high  frequency  electron  plasma  wave,  (ol  is  the  so-called  ion 

sound  frequency  and  ncr  is  the  critical  density  in  the  plasma  at  which  an  incident  light 
wave  is  reflected. 

3.    Self  Generated  Magnetic  Fields  in  Plasmas 

Experiments  and  calculations  have  shown  the  existence  of  self-generated 
magnetic  fields  [Ref.  3].  These  magnetic  fields  can  be  generated  in  many  different  ways 
[Ref.  4].  Let's  consider  the  simplest  way.  Consider  a  target  irradiated  by  high  intensity 
laser  light.  There  will  be  both  a  density  and  temperature  gradient  in  the  resulting  plasma, 
as  indicated  in  Figure  1 .  When  the  cross  product  between  the  density  gradient  ( Vn)  and 
the  temperature  gradient  (VT)  is  nonzero,  a  magnetic  field  is  generated.  Faraday's  law 

shows  that  B  is  nonzero  (i.e.,  a  magnetic  field  grows  up)  whenever  there  is  a  curl  of  the 
electric  field.  In  an  expanding  plasma,  the  electric  field  is  proportional  to  the  gradient  of 
the  pressure  (VnKT)  divided  by  the  density.  One  then  finds  that  a  magnetic  field  is 
generated  whenever  Vn  x  VT  *  0  . 


Target 


VT 


Vn 


Figure  1.  A  schematic  illustrating  crossed  density  and  temperature  gradients 
generated  when  a  target  is  irradiated. 


B.  COMPUTER  SIMULATION  OF  PLASMA  USING  PARTICLE  CODES 

A  theoretical  description  is  often  used  where  the  plasma  is  treated  as  two  charged 
fluids.  We  will  use  a  complementary  numerical  description  of  plasma  behavior  using 
particle  codes  [Ref.  1,5,6].  Computer  simulation  of  plasma  using  particle  codes  is  a  very 
direct  and  powerful  approach,  particularly  for  investigating  kinetic  and/or  nonlinear 
effects.  Since  we  are  investigating  strongly  nonlinear  plasma  effects,  we  utilize  the 
power  of  the  computer  to  simulate  the  behavior  of  the  charged  particles  of  plasma. 
The  approach  is  as  follows:  numerically  follow  the  motion  of  a  large  collection  of 
charges  in  their  self-consistent  electric  and  magnetic  fields.    The  basic  cycle  is  illustrated 
in  Figure  (2). 

From  the  positions  and  velocities  at  any  given  time,  we  compute  the  charge  and 
current  densities  on  a  spatial  grid  that  is  fine  enough  to  resolve  the  collective  behavior  of 
the  charges.  Using  the  computed  charge  and  current  densities,  then  compute  the  self- 
consistent  electric  and  magnetic  fields  via  Maxwell's  equations.    Next,  we  use  these 


fields  in  the  equations  of  motion  to  advance  the  positions  and  velocities  of  the  charges. 
Finally,  continue  around  this  basic  cycle  with  a  time  step  sufficiently  small  to  resolve  the 
highest  frequency  in  the  problem  (which  is  often  the  electron  plasma  frequency  (ope). 


Newton's  /  \        Grid 

laws  /  \       assignment 


{£•6}  {f.j} 


Maxwell's 
equations 

Figure  2.  The  Basic  Cycle  of  a  Particle  Simulation  Code,  From  Ref.[l]. 

While  plasma  actually  consists  of  a  huge  number  of  particles,  particle  simulations 
are  limited  by  the  capabilities  of  computers  to  104  to  108  particles,  and  a  single  simulation 
particle  must  therefore  represent  a  large  number  of  actual  particles.  Since  we  seek  to 
describe  the  collective  behavior  of  the  plasma,  we  can  follow  behavior  on  characteristic 
length  £ ,  which  is  larger  than  the  Debye  length  XDe ,  to  simplify  our  analysis  and  not 
consider  the  fine-scale  fluctuating  microfields  associated  with  particle  discreteness.  This 
simplification  allows  us  to  represent  the  fields  on  a  discrete  grid  in  space,  which  is  fine 
enough  to  resolve  the  collective  behavior  and  sufficiently  coarse  to  ignore  the 
microfields.  A  mesh  size  8  «  ADe  is  usually  chosen. 

Alternatively,  this  spatial  averaging  process  can  be  viewed  as  using  "finite-size" 
particles  of  size  5  «  XDe .  The  behavior  of  a  group  of  such  particles  is  the  same  (with 


minor  modifications)  as  the  behavior  of  point  charges  for  scale  lengths  I »  8,  but 
fluctuations  with  scale  lengths  £  «  8  are  suppressed.  Hence,  we  achieve  with  a  "trick" 
what  nature  does  with  the  use  of  an  enormous  number  of  particles;  i.e.,  smooth  over  the 
fine  scale  length  microfields. 

1.    Basic  Ingredients  of  a  Particle  Code 

All  particle  codes  share  some  basic  ingredients,  which  we  shall  illustrate.  We 
treat  the  ions  as  a  fixed,  neutralizing  background,  assume  no  imposed  magnetic  fields, 
and  consider  only  electrostatic  fields.  In  this  electrostatic  limit,  the  magnetic  field 
generated  by  plasma  currents  is  negligible,  and  Maxwell's  equations  reduce  to  Poisson's 
equation,  which  for  simplicity  we  look  at  in  only  one  dimension: 

f0VV  =  £0  df/dx2  =  -e{nt  -  ne)  (2) 

In  this  equation,  So  is  the  permittivity  constant  of  free  space,  <f>  is  the  electric  field 

potential,  e  is  the  electron  charge,  ne  and  nl  being  the  electron  and  ion  charge  densities 

respectively.  Variations  are  allowed  in  only  one  direction,  and  periodic  boundary 

conditions  are  adopted. 

Our  first  task  is  to  compute  the  charge  density.  Our  system  extends  from  0  to  L, 

as  shown  in  Figure  (3).  We  divide  this  system  into  NC  cells  and  for  convenience  take  the 

cell  size  (8)  to  be  unity  (i.e.,  L  =  NC).  The  grid  points  are  identified  as  the  integer  values 

of  the  position,  augmented  by  one  so  that  the  counting  begins  with  one.  Note  that  the 

(NC  +  1)*  grid  point  is  then  identical  with  the  first,  due  to  the  assumed  periodic  boundary 

conditions.  Given  the  position  of  the  charge,  many  schemes  can  be  used  to  assign  it  to 

the  spatial  grid.  For  example,  we  could  just  assign  the  charge  to  its  nearest  grid  point 


location.  A  better  scheme  is  to  share  the  charge  between  its  two  nearest  grid  points,  (see 
Figure  3) 


x=0 


Particle  i 
/  J 


x=L 


ceil  1  celli  +  l 

Figure  3.  The  charge-sharing  scheme  for  finite  sized  particles,  From  Ref.[l]. 


For  a  charge  located  a  distance  A  x  to  the  right  of  the  i*  grid  point,  we  then  have: 

Ap(i)  =  q(l-Ax)  (3) 

Ap(i+l)  =  qAx  , 
where  Ap  is  the  increment  to  the  charge  density. 

Having  assigned  the  charges  and  determined  the  charge  density  on  the  spatial 
grid,  we  next  determine  the  self-consistent  electric  field  using  Poisson's  equation 

dE/dx  =  p/e0  (4) 

The  simplest  approach  is  to  finite-difference  Poisson's  equation, 

E{i  +  1)  =  E(i)  +  S[(p{i  +  1)+  p(i))/2l\/e0)  (5) 


An  alternative  approach  is  to  Fourier  transform  the  charge  density,  use  Poisson's 
equation  to  find  Ek  (ikEk  =  pk  js) ,  and  invert  the  transform  to  determine  the  electric  field 
on  the  spatial  grid. 

The  last  step  in  the  basic  cycle  is  to  use  the  electric  field  to  move  the  particles. 
The  electrical  force  is  assigned  from  the  grid  to  the  individual  particles  using  the  same 
scheme  chosen  to  assign  the  charges  to  the  grid.  For  example,  considering  a  particle  a 
distance  Ax  to  the  right  of  the  Ith  grid  point  and  using  linear  interpolation,  we  obtain  for 
the  force  F  on  the  particle 

F  =  ?£(0[l  -  Ax]+ qE{i  +  l)Ax  (6) 

In  this  code,  Poisson's  equations  are  only  solved  for  the  electric  field  and  the  magnetic 
field  is  assumed  zero.  Therefore,  the  velocity  (v)  and  position  (x)  of  each  particle  can  be 
advanced  At  in  time  using  a  "leap  frog"  algorithm;  i.e., 

v"+1/2=v"-1/2+F"Ar  (7) 

x"+l  =xn+v"+l/2At  (8) 

The  superscripts  denote  time  step.  By  defining  x  and  v  one-half  time  step  apart,  we 
achieve  second-order  accuracy  in  the  time  step.  In  the  initial  conditions,  x  and  v  are 
defined  at  the  same  time  (t=0),  but  it  is  straightforward  to  then  displace  the  velocities 
backward  in  time  using  the  force  at  t=0. 

The  plasma  evolution  is  computed  by  simply  continuing  around  the  basic  cycle 
using  a  time  step  small  enough  to  resolve  the  characteristic  oscillations  of  the  plasma.  In 
this  electrostatic  limit,  the  highest  frequency  oscillation  has  a  frequency  near  cope,  the 
electron  plasma  frequency.  Hence,  a  time  step  of  about  0.2o)~l  is  commonly  used 
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Of  course,  it  is  also  necessary  to  resolve  the  scale  lengths  characteristic  of  the  collective 
behavior.  Hence  a  grid  size  of  about  X,de  is  commonly  used.  It  is  interesting  to  note  that  a 
numerical  instability  would  be  introduced  if  the  grid  size  were  chosen  to  be  too  many  A,de. 
This  instability  is  due  to  aliasing,  which  arises  from  the  fact  that  on  the  grid  a  disturbance 
with  wave  number  k  cannot  be  distinguished  from  spurious  ones  with  wave  numbers  k  + 
27m/ S  where  8  is  the  grid  spacing  and  n  is  an  integer  [Ref.  1,19]. 

2.  PIC  Code  Programming  Considerations 

Random  particle  initialization  is  the  simplest  and  most  commonly  used  method. 
However,  it  introduces  noise  in  the  computations  which  can  mask  interesting  physical 
phenomena,  particularly  when  we  wish  to  study  the  low  amplitude  linear  regime.  In  such 
cases  the  programs  utilize  what  is  referred  to  as  a  "quiet  start"  which  avoids  this 
background  noise  introduced  by  the  randomness  of  the  numbers.  The  program  that 
simulates  the  plasma  particle  behavior  can  remain  in  this  "quiet"  state  for  only  a  finite 
time,  which  must  be  long  enough  to  allow  observation  of  the  low-amplitude  phenomena 
being  studied.  In  such  cases,  staggering  the  particles,  so  that  they  do  not  form  exact 
beams,  is  an  effective  way  to  delay  the  onset  of  the  noisy  behavior. 

Even  with  this  simplest  of  particle  code,  many  instructive  and  interesting  plasma 
problems  can  be  examined.  As  one  example,  let's  use  ESI,  a  flexible  and  optimized  Id 
electrostatic  PIC  code  to  look  at  the  strongly  non-linear  behavior  of  a  two-stream 
instability  [Ref.  5]. 
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In  this  simple  case,  half  the  electrons  are  given  an  initial  drift  velocity  in  one 
direction  and  half  in  the  other  direction.  This  counter  streaming  motion  is  a  source  of 
free  energy,  which  drives  unstable  a  variation  in  density  and  potential. 

Figure  4a  shows  the  electron  phase  space  (velocity  versus  position  for  all  the 
electrons),  electric  field  and  the  potential  at  an  early  time.  Note  the  counter  streaming 
plasmas.  In  Figures  4b  and  Figure  4c,  the  effect  of  the  instability  on  the  electric  field  and 
the  potential  in  the  plasma  can  be  seen.    Finally,  the  growing  electric  field  becomes  so 
large  that  vortices  form  in  phase  space,  as  the  electrons  are  alternatively  slowed  down 
and  speeded  by  the  large  electric  field,  as  shown  in  Figure  4c.    Figure  4d  shows  time 
histories  of  the  kinetic  energy,  the  electric  field  energy  and  the  total  energy  of  the  system. 
Note  the  exponential  growth  of  the  electric  field  energy  at  early  times  and  the  subsequent 
saturation.  As  the  electric  field  energy  increases,  the  kinetic  energy  decreases.  The  total 
energy  varies  very  little  compared  to  these  changes,  showing  that  the  code  is  conserving 
energy  with  sufficient  accuracy. 
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Figures  4a-c  show  Vx  Phase  Space  ,  Electric  Field,  and  Potential  respectively 
versus  x.  Figure  4  d  shows  Kinetic,  Electric  Field,  and  Total  energies  respectively  versus 
time. 

Zohar  represents  a  significant  improvement  to  ESI  in  PIC  modeling  of  plasma 
behavior.  This  code  solves  all  of  Maxwell's  equations  in  order  to  include  light  waves  and 
uses  equations  of  motion  that  are  relativistically  correct.  Zohar  (which  means  brilliance 
or  radiance  in  Hebrew)  is  the  2  14  D  PIC  code  used  at  Lawrence  Livermore  Labs  to  model 
the  interaction  of  intense  laser  light  with  plasmas  and  for  some  electron  beam-plasma 
studies[Ref.  7].  2  XA  D  means  that  variations  are  followed  in  two  directions,  and  that  all 
three  velocities  are  included  in  the  code  calculations.  We  note  that  state  of  the  art,  three- 
dimensional  PIC  codes,  which  run  on  parallel  machines,  now  exist  [Ref.  8,9].  The 
simulations  that  follow  are  results  of  Zohar  runs. 
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II.  ULTRA-INTENSE  REGIME  OF  LASER  MATTER  INTERACTIONS 

The  recent  development  of  Ultra  Short  Ultra-Intense  (US-UI)  laser  pulses  has 
allowed  for  the  generation  of  laser  intensities  that  are  orders  of  magnitude  over  what  was 
possible  with  previous  laser  technology  [Ref.10].  The  ultrashort  regime  includes  laser 
light  with  pulse  lengths  xp  of  about  10  femtoseconds  (fs)  up  to  10  picoseconds  (ps).  The 
ultra-intense  regime  covers  intensities  starting  at  about  1017  W/cm  .  The  ultra-intense 
regime  is  usually  defined  as  the  point  at  which  the  normalized  amplitude  of  the  electric 
field  of  the  laser  light  is  greater  than  or  equal  to  1,  i.e.,  eEo/me0)dC  >  1  [Ref.l  1].  Here, 
E0is  the  electric  field  amplitude  of  the  laser  light,  me  is  electron  rest  mass,  co0  is  the 
frequency  of  the  laser  light,  and  c  is  light  speed.  Even  as  late  as  ten  years  ago,  most 
lasers  could  not  produce  light  that  could  approach  this  electric  field  strength,  with  the 
exception  of  a  few  CO2  lasers  built  in  the  1970's.  Now,  there  are  a  large  number  of  ultra- 
intense  lasers  all  over  the  world.  What  made  this  possible  was  the  application  of  chirped 
pulse  amplification  (CPA)  [Ref.  10],  originally  developed  in  microwave  technology,  to 
laser  systems.  This  breakthrough  in  laser  hardware  has  spawned  considerable  interest  in 
US-UI  laser-matter  interactions.  Never  before  has  it  been  possible  to  deposit  so  much 
laser  energy  in  such  a  short  amount  of  time,  in  such  a  tiny  volume  on  the  front  of  a  solid 
target. 

Chirped  Pulse  Amplification  is  critical  because  laser  pulses  of  extremely  high 
power  density  (gigawatts/centimeter2)  can  severely  damage  optical  components  such  as 
amplifiers,  lenses,  and  mirrors.  CPA  technology  allows  very  short  laser  pulses  to 
generate  extremely  high  peak  powers  by  stretching  the  initially  low-energy  laser  pulse 
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more  than  10,000  times  its  duration  before  amplifying  it  and  recompressing  the  pulse 
back  to  near  the  original  duration  after  amplification.  Damage  to  the  laser  optics  is 
avoided  by  passing  the  long  low  power  pulse  through  the  amplifier  (see  Figure  5). 


lni£a'  sho-i  pulse 


Sftei-pulss  osa";alo- 


Power  ampJifie.-s  /fa      «■      •  ■  »       /£_ w 


Hich-ensfgy 
Otirashori  pulse 


Figure  5.    The  concept  of  chirped-pulse  amplification  and  other  new  technologies 
enabling  the  production  of  ultra-intense  pulses,  From  Ref.[10]. 

Previously,  the  intense  CO2  lasers  were  nanoseconds  long,  which  creates  large 
preformed  plasmas  that  the  laser  interacted  with.  There  was  then  a  considerable  amount 
of  absorption  in  the  underdense  part  of  the  plasma,  as  well  as  enough  time  for  the  laser  to 
filament  on  its  way  to  the  critical  surface.  This,  in  turn,  led  to  a  large  fraction  (>  10%)  of 
the  laser  energy  converted  to  hot  electrons,  which  were  thought  to  be  created  by  various 
laser-plasma  interactions.  However,  US-UI  laser  pulses  can  have  pulse  lengths  of  the 
order  of  picoseconds  and  shorter.  This  means  that  they  have  the  ability  to  both  create  a 
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dense  plasma  (sometimes  several  times  the  critical  density  and  up  to  solid  densities)  and 
interact  with  it,  if  proper  care  is  taken  to  eliminate  the  pre-pulse. 

The  interaction  of  an  intense  laser  with  very  dense  plasma  is  interesting  for 
several  reasons.  First,  since  the  laser  light  is  absorbed  most  efficiently  near  the  critical 
surface,  the  material  can  potentially  absorb  a  large  fraction  of  the  laser  light.  Laser  light 
transfers  energy  to  the  plasma  via  the  electrons.  Electrons  readily  oscillate  in  the  electric 
field  of  the  laser.  This  motion  is  referred  to  as  the  quiver  velocity,  and  is  defined  as 

P0Jmec  =  rv0Jc=  eE0/mecco0  =  Ju2jl3x\0lt  (9) 

Here,  posc  (vosc)  is  the  transverse  quiver  momentum  (velocity)  of  an  electron  in  the  laser 
field  with  an  electric  field  amplitude  of  Eo,  I  the  laser  intensity  (in  W/cm2),  X^  is  the  laser 
wavelength  in  microns,  and  y  is  the  relativistic  factor  defined  below.  The  energy 
associated  with  this  momentum  can  easily  reach  MeV  levels  for  electrons  in  the  presence 
of  lasers  with  (normalized)  intensities  101  Wum2/cm2.  Therefore,  relativistic  effects 
become  important  for  these  interactions,  and  are  included  in  calculations  concerning  US- 
UI laser-matter  interactions.  Note  that  the  quiver  velocity  scales  in  terms  of  IX.2.  This  is 
a  common  scaling,  and  so  it  is  useful  to  introduce  this  as  the  normalized  intensity.  In  this 
way,  many  results  discussed  can  easily  be  applied  to  any  wavelength  laser  desired. 
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A.  SIMPLE  MODELS  FOR  REALISTIC  LASERS  AND  SOLIDS 

An  ideal  model  for  ultra-intense  pulses  interacting  with  solids  would  require  that  a 
perfect  Gaussian  laser  pulse,  with  a  temporal  pulse  length  (FWHM)  of  t0  and  a  peak 
intensity  of  Iq,  strike  a  solid  with  a  step-like  density  interface  [Ref.l  1].  When  the  laser 
intensity  exceeds  about  1012  W/cm2,  ionization  of  the  electrons  in  the  target  occurs,  and  a 
plasma  is  formed  on  the  surface  of  the  solid.  Typical  ion  densities  are  n;  »  1023  cm"3, 
and  electron  densities  are  rie  =  Z  n;.  For  US-UI  pulses,  this  happens  quite  early  in  the 
interaction,  and  so  most  of  the  pulse  interacts  with  a  plasma,  as  opposed  to  an  un-ionized 
solid.  The  relativistically  correct  dispersion  relation  (for  circular  polarization)  that 
governs  the  propagation  of  light  waves  in  the  plasma  is  given  by 


kc  =  4co2pely-G>l  (10) 

where  k(co0)  is  the  wavenumber  (frequency)  of  the  laser  light, 

r  =  Ji  +  P2JM  (ii) 

is  the  relativistic  factor,  and 

cope  =  ^4mee2/me  (12) 

is  the  plasma  frequency. 

First,  consider  the  nonrelativistic  case  where  y  «  1 ,  which  is  valid  up  to 
(normalized)  intensities  of  about  5  x  1017  W-um2/cm2.  In  this  case,  k  approaches  zero 
when  G)pe=  ycD0 .  The  laser  no  longer  propagates,  but  is  reflected  when  this  condition  is 
satisfied.  The  density  at  which  this  occurs  is  called  the  critical  density  and  is  given  by 

ncr=yl.U1021/AMcm-3  (13) 
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where  X^  is  the  laser  wavelength  in  micrometers.  The  laser  light  does  penetrate  slightly 
past  critical,  but  falls  off  exponentially  in  a  distance  determined  by  the  plasma  density. 
This  distance  is  known  as  the  collisionless  skin  depth  and  is  given  by  c/o)pe  in  the 
nonrelativistic  limit.  For  example,  for  a  solid  target  of  carbon  (Z  =  6),  the  skin  depth  is 
about  6  nm  (nanometer  =  10"9  m). 

For  the  relativistic  case  (y  >1),  we  see  that  the  effective  density  at  which  the  light 
will  be  reflected  from  the  plasma  has  increased  by  the  factor  y.  Thus,  in  the  relativistic 
case  laser  light  can  propagate  even  further  into  the  overdense  plasma  due  to  the  mass  of 
the  electrons  increasing  as  they  become  relativistic.  This  increase  in  inertia  makes  it 
more  difficult  to  generate  the  oscillating  current  in  the  plasma  needed  to  reflect  the  laser 
light. 

A  second  interesting  effect  evident  in  US-UI  laser  matter  interactions  is  that  the 
laser  pressure  can  easily  be  larger  than  the  plasma  pressure,  even  for  extremely  dense 
plasmas.  The  pressure  associated  with  a  laser  with  intensity  I  (in  W  /cm2)  is 

pL  =3307/1 018Mbar  (14) 

Efficient  coupling  of  laser  light  to  electrons  has  been  exploited  in  two  ways  to 
date.  First,  there  has  been  considerable  effort  at  trying  to  use  these  interactions  to 
produce  high  brightness  X-ray  lasers  [Ref.  11].  Second,  the  possibility  of  generating 
large  numbers  of  energetic  electrons  has  been  proposed  as  a  novel  way  to  deliver  a  large 
amount  of  energy  to  a  hot  spot  in  a  compressed  core  of  deuterium  and  tritium  in  order  to 
ignite  a  fusion  pellet.  This  is  commonly  known  as  the  Fast  Igniter  fusion  scheme 
[Ref.  12].  Central  to  both  of  these  applications  is  the  question  of  efficiency.  How 
efficient  is  the  coupling  of  laser  light  to  the  production  of  hot  electrons?  How  does  it 
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scale  with  laser  intensity,  wavelength,  f-number,  or  spot  size?  How  does  the  efficiency 
scale  with  the  plasma  parameters  such  as  temperature,  density,  and  scale  length?  [Ref.13] 
Computer  simulations  using  Particle-In-Cell  codes  were  used  to  provide  an  early 
glimpse  into  the  ultra-intense  regime.  As  a  simple  example,  consider  a  two-dimensional 
simulation  of  a  focused  light  beam  with  intensity  IAM  =  1 019  W  ■  fjm2 1  cm2  incident  onto 

an  overdense  plasma.  Over  a  third  of  the  light  was  found  to  be  absorbed  into  very 
energetic  electrons  as  shown  by  Figure  6a  [Ref.  13].    These  electrons  had  an  effective 
temperature  of  about  1  MeV  at  this  intensity.  As  shown  in  Figure  6b,  the  light  pressure 
punched  a  hole  in  the  plasma,  causing  it  to  recede  at  a  significant  fraction  of  the  velocity 
of  light.  Both  of  these  effects  have  now  been  observed.  Very  large  self-generated 
magnetic  fields  were  found  in  these  simulations  (up  to  109  Gauss!). 
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Figure  6a,  From  Ref.[13].  A  plot  of  the  energy  distribution  of  escaping  electrons. 
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Figure  6b,  From  Ref.[13].  A  plot  of  electrons  in  space  showing  a  hole  being  punched 
into  the  plasma. 
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III.    EFFECT  OF  LARGE  SELF-GENERATED  MAGNETIC  FIELDS  ON 
ULTRA  INTENSE  LASER  PLASMA  INTERACTIONS 


The  focus  of  this  thesis  is  to  show  how  large  self-generated  magnetic  fields  affect 
the  interaction  of  ultra-intense  laser  pulses  with  plasmas  as  found  in  the  early  simulation 
studies  just  discussed,  the  magnitude  of  the  self-generated  magnetic  (BdC)  field  can  be 
estimated  as  Bos/3,  where  Bos  is  the  amplitude  of  the  oscillating  magnetic  field  of  the 
short,  intense  laser  pulse.  In  other  words,  about  10%  of  the  energy  of  the  oscillating 
magnetic  field  was  "rectified".  When  targets  are  irradiated  with  1.06um  laser  light  with 
intensities  ranging  from  1018  W/cm2  to  1021  W/cm2,  Bdc  can  become  very  large,  with 
typical  values  of  30-900  MG  as  shown  in  Table  I.  To  put  such  large  magnetic  fields  into 
perspective,  we  note  the  following: 

-  The  approximate  magnetic  field  at  the  earths  surface 

*  10-4 T*l  Gauss. 

-  The  approximate  magnetic  field  near  a  big  electromagnet 

*  1.5  T*  1.5  xlO4  Gauss. 

-  The  approximate  magnetic  field  near  the  surface  of  a  neutron  star 
«108T*1012  Gauss  [Ref.14]. 
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I(W/cm2) 

Bdc(MG) 

7 

CDce/CDo 

101* 

30 

1.2 

0.7 

10,y 

90 

2.1 

0.4 

1<T 

300 

6.1 

0.5 

10*' 

900 

19 

0.5 

Table  1.  Typical  estimated  values  of  Bdc  and  cyclotron  frequency  in  the  ultra- 
intense  laser  regime.  Here  ooce  =  eBdc/mcY  and  y  =  (l  +  IA2M  /2.8  x  1 018 ) 

Some  other  characteristic  parameters  are  also  shown  in  Table  I.  The  electron 
motion  in  such  intense  laser  light  can  become  relativistic.  The  relativistic  factor  y  is 
estimated  as: 

7  =  (l  +  /4/2.8xl018}^,  (15) 

where  linearly  polarized  light  has  been  taken.  Typical  values  of  y  range  from  1 .2  to  19. 
The  characteristic  frequency  with  which  electrons  gyrate  (see  Figure  7)  about  the 
magnetic  field  is  the  so-called  electron  cyclotron  frequency  coce,  where 

<?ce=eBdc/my  (16) 

As  shown,  CDce  becomes  a  significant  fraction  of  coo  ,  the  laser  light  frequency. 
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Figure  7.  Cyclotron  motion  of  an  electron. 

Large  self-generated  fields  can  significantly  modify  the  interaction  of  the  laser 
light  with  the  plasma  [Ref.  15].  This  phenomenon  can  be  explained  by  considering  the 
oscillatory  motion  of  the  electrons,  which  results  from  the  electric  field  of  the  incident 
laser  light.  The  Lorentz  force  associated  with  this  oscillation  of  electrons  across  the 
magnetic  field  is  determined  by  : 

F  =  -euosxBdc  .  (17) 

In  this  equation,  uos  is  the  velocity  of  electron  oscillation  in  the  laser  field.  This  Lorentz 
force  can  become  very  strong.  For  Bdc  =  109  Gauss  and  taking  the  peak  uos  =  c,  the  force 
corresponds  to  an  oscillating  electric  field  with  amplitude  of  approximately  3xlOu  V/cm! 
This  large  force  can  be  expected  to  enhance  the  absorption  efficiency,  increase  the  heated 
electron  energies,  and  lead  to  the  generation  of  harmonics  of  the  incident  light.  For  a 
discussion  of  some  of  these  effects  at  much  lower  intensities,  see  Ref.  16. 

To  demonstrate  the  effects  of  the  self-generated  fields,  we  used  so  called 
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1  1/2  D  simulations  (1  position,  2  velocities  i.e.  velocities  in  two  directions  for  electrons 
and  ions)  with  Zohar,  a  Particle-In-Cell  code  (PIC)  that  solves  the  complete  set 
Maxwell's  equations  and  uses  relativistic  dynamics  [Ref.  7]. 

As  shown  in  Figure  14,  an  intense  light  pulse  was  sent  from  vacuum  onto  a  slab  of 
overdense  plasma. 
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Figure  14.  A  schematic  illustration  of  a  computer  simulation  of  intense  laser  light 
interacting  with  an  overdense  plasma. 


In  these  simulations,  a  dc  magnetic  field  pointing  in  the  y  direction  was  added  in  the 
code.  The  oscillating  electric  field  of  the  laser  light  Eos  was  chosen  to  be  either  parallel 
or  perpendicular  to  Bdc  (  by  simply  rotating  the  electric  vector  of  the  laser  light  90°  ). 
Some  typical  parameters  are  shown  in  this  figure.  The  initial  plasma  density  rose  rapidly 
to  50x  ricr  where  ricr  is  the  critical  density.  The  intensity  of  incident  light  rose  to  its  peak 
and  then  returned  to  zero  in  a  time  of  360  co^1 .  This  corresponds  to  a  pulse  width  of  1 80 
fs  (femtoseconds  or  10"15  seconds)  for  1.06um  light. 
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To  most  simply  illustrate  the  effects,  consider  an  example  in  which  the  1 .06  urn 
laser  light  intensity  is  1020  W/cm2,  the  magnitude  of  the  imposed  Bdc  is  300  MG,  and  the 
ions  are  fixed.  The  enhanced  heating  due  to  the  dc  magnetic  field  is  apparent  in  Figure 
(15),  which  shows  electron  phase  space  from  simulations  with  and  without  the  external 
Bdc  field.  More  electrons  are  heated  to  higher  velocity  with  the  B  field,  as  expected  due 
to  the  additional  force  (-e  uos  x  BdC).  Similar  results  can  be  seen  from  the  distribution 
functions  (number  of  electrons  vs.  energy)  shown  in  Figure  (16). 
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Figures  15  a,b.  Snapshots  of  electron  phase  space  in  simulations  a)  with  Bdc  and  b) 
without  Bdc- 
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Figures  16  a,  b.  Plots  showing  numbers  of  electrons  versus  energy  a)  with  Bdc  field 
on  b)  without  BdC  field. 


Absorption  can  be  measured  by  several  different  methods.  In  our  approach,  we 
chose  to  measure  absorption  by  computing  the  change  in  kinetic  energy  of  the  plasma 
(from  times  before  and  after  the  laser  pulse)  and  then  dividing  by  the  energy  of  the  laser 
pulse.  Figure  1 7  shows  the  absorbed  energy  versus  time.  In  these  simulations,  the  upper 
curve  corresponds  to  the  simulation  with  the  magnetic  field  and  the  lower  curve  is  the 
simulation  without  the  magnetic  field.  If  we  normalize  to  the  energy  of  the  laser  pulse, 
the  fractional  absorption  was  approximately  19%  without  Bdc;  it  rose  to  about  30%  with 
Bdc- 
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We  also  ran  a  simulation  with  BdC,  but  with  it  pointing  in  the  Z-direction.  In  this 
orientation  uos  and  Bdc  are  parallel,  meaning  that  u0s  x  Bdc  =  0.  As  we  expected,  there  was 
no  enhancement  of  the  absorption  observed. 
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Figure  17  showing  enhanced  absorption  when  BdC  is  turned  on. 

The  simulations  also  show  that  the  force  due  to  oscillations  of  electrons  across  Bdc 
leads  to  an  intense  generation  of  even  harmonics  in  the  reflected  light  (see  Figure  18). 
These  even  harmonics  arise  from  the  beating  of  the  incident  light  wave  with  the  density 
fluctuation  driven  by  the  uos  x  B^  force.  In  addition,  there  is  enhanced  radiation  at  lower 
frequency  i.e.  co  <  coce  which  we  attribute  to  radiation  by  electrons  as  they  gyrate  about 

the  strong  Bdc  field.  It  would  be  interesting  to  look  for  this  radiation  in  experiments,  as 
detection  of  this  radiation  may  serve  as  a  diagnostic  for  the  large  self-generated  fields  that 
we  propose  can  be  generated  in  the  ultra-intense  laser  regime. 
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Figure  18.  Frequency  components  of  the  reflected  electric  field  showing  even 
harmonics  with  Bdc. 

Lastly,  we  carry  out  the  same  simulations  allowing  the  ions  to  move.  The 
electron-ion  mass  ratio  is  taken  to  be  1/3700  (i.e.  nVmi  =  1/3700).  As  shown  in  Figure 
19,  we  again  find  that  the  absorption  is  greater  with  Bdc  («  20%)  than  without  Bdc  (fabs~ 
2%).  It  is  somewhat  reduced  in  both  cases  because  the  intense  pressure  of  the  incident 
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light  wave  strongly  steepend  the  density  profile  of  the  plasma.  Again,  we  find  that  the 
heated  electrons  are  significantly  more  energetic  in  the  simulation  with  Bdc.  This  is  clear 
from  the  electron  phase  spaces  shown  in  Figure  20. 
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Figure  19.  Plots  showing  enhanced  absorption  with  Bdc. 
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Figure  20.  Plots  of  electron  phase  space  in  simulations  with  moving  ions  a)  with 
and  b)  without  Bdc. 
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The  ion  phase  space  shown  in  Figure  21  illustrates  the  ion  acceleration  into 
vacuum  on  both  sides  of  the  plasma.  This  acceleration  is  produced  by  the  electric  field 
that  sets  up  to  prevent  the  electrons  from  leaving  the  ions  behind.  Note  also  in  Figure  21 
the  strong  ion  acceleration  as  the  pressure  of  the  intense  light  pushes  on  the  plasma.  With 
Bdc,  ions  that  are  more  energetic  are  generated  both  in  the  expanding  plasma  and  at  the 
light-plasma  interface. 
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Figure  21.  Simulations  of  Ion  phase  space  with  and  without  Bdc 
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Finally,  the  simulations  show  a  very  interesting  compression  of  the  Bdc.  The 
pressure  of  the  intense  light  wave  compresses  the  density  of  the  plasma,  as  it  reflects  (and 
partially  absorbs).  The  peak  magnitude  of  the  field  can  be  significantly  increased.  As 
seen  in  Figure  22,  the  peak  Bdc  has  gone  up  by  a  factor  of  about  three,  which  corresponds 
to  over  109  Gauss  in  this  particular  example. 
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Figure  22.  Simulation  of  the  spatial  profile  of  the  magnitude  of  the  self-generated 
magnetic  field  after  irradiation  by  an  intense  laser  pulse  with  intensity  of  1021 
W/cm2. 
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IV.  CONCLUSION 

We  have  investigated  via  computer  simulations  how  self-generated  magnetic 
fields  up  to  109  Gauss  affect  the  interaction  of  short,  ultra-intense  pulses  of  laser  light 
with  plasmas.  We  have  shown  that  the  laser  absorption  and  the  high-energy  electron 
generation  can  be  significantly  enhanced  by  these  self-generated  magnetic  fields.  The 
enhancement  is  attributed  to  the  very  large  Lorentz  forces  resulting  from  the  oscillation 
of  electrons  by  the  laser  light  across  these  fields. 

Other  effects  shown  are  the  generation  of  intense  harmonics  in  the  reflected  light 
and  noticeable  cyclotron  radiation.  Finally,  the  large  self-generated  magnetic  fields  can 
be  compressed  to  even  large  values  by  the  pressure  of  the  ultra-intense  laser  light. 

Since  the  late  1970's,  the  existence  of  self-generated  magnetic  fields,  in  laser- 
irradiated  targets,  have  been  investigated.    Some  of  the  initial  research  confirming  the 
existence  of  these  magnetic  fields  was  conducted  by  scientists  at  the  Naval  Research  Lab 
[Ref.  17].  To  attain  magnetic  fields  on  the  order  of  109  Gauss,  which  we  have 
considered,  requires  a  powerful  laser  capable  of  irradiating  a  target  with  intensities  on 
the  order  of  1021  W/cm2.     Already,  such  intensities  have  been  achieved  by  the 
"Petawatt",  short  pulse  laser  at  the  Lawrence  Livermore  National  Laboratory. 

Exciting  high-energy  applications  and  new  research  in  creating  Antimatter  using 
ultra-intense  lasers  like  Petawatt  laser  is  being  conducted.  As  an  example,  intense  light 
from  the  Petawatt  laser  has  been  directed  onto  a  thin  gold  film  where  it  creates  a  plasma 
plume,  which  acts  as  a  sort  of  messy  Wakefield  accelerator  [Ref.  18]  (see  Figure  23).  In 
particular,  the  electric  field  of  the  laser  light  rips  electrons  from  the  gold  atoms  and  send 
the  electrons  shooting  off  with  energies  as  high  as  100  MeV. 
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Figure  23.  Illustration  of  an  ultra-intense  laser  matter  interaction.  The  laser's 
electric  field  is  strong  enough  to  rip  electrons  from  the  gold  atoms  and  shoot  them 
off  with  energies  up  to  100  MeV,  initiating  the  generation  of  y  rays,  which  in  turn 
can  initiate  fission  or  the  generation  of  anti-matter  (electron-positron  pairs),  From 
Ref.[18  ]. 
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The  generation  of  the  relatively  lower  3-10  MeV  electrons  using  ultra-intense 
laser  light  (intensity  about  1021  W/cm2),  is  an  important  issue  for  laser  radiography. 
Some  of  the  higher  energy  electrons  radiate  gamma  rays,  which  in  turn  can  induce  fission 
or  create  electron-positron  pairs  (the  first  antimatter  made  in  laser-solid  interactions).  A 
positron  is  an  electron  with  a  positive  charge  and  is  the  anti-particle  of  the  electron.  The 
positron  does  not  exist  naturally  on  our  planet.  It  'arrives'  in  cosmic  radiation  or  is 
created  in  the  laboratory.  Thus,  laser  photons  at  the  electron-volt  level  can,  by  teaming 
up,  initiate  the  sort  of  million-electron-volt  nuclear  reactions  that  normally  take  place 
only  in  an  accelerator.  Moreover,  the  femtoseconds  laser  pulses  can  be  focused  to  a 
much  smaller  spot  size  than  is  possible  with  any  conventional  particle  beam  . 

Continued  enhancements  and  improvements  in  technology  offer  rich  and  exciting 
frontiers  of  science  to  be  explored.  Continued  probing  and  exploring  of  the  ultra-intense 
regime  of  laser  plasma  interactions  shall  give  more  understanding  of  very  high 
temperature  matter.  This  understanding  will  define  and  point  to  new  applications, 
ranging  from  developments  in  laser  radiography  to  the  successful  application  of  the  fast 
ignition  concept  of  inertial  confinement  fusion. 


37 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


38 


LIST  OF  REFERENCES 

1 .  Kruer,  William  L.,  Physics  of  Laser  Plasma  Interactions,  (Addison- Wesley), 
Redwood  City,  1988 

2.  Chen,  Francis  F.,  Introduction  to  Plasma  Physics  and  Controlled  Fusion  Vol  I. 
(Plenum  Press),  New  York,  NY,  1984 

3 .  Stranger,  J.  A.  et  al.,  Phys.  Rev.  Letters  26,  1 0 1 2  ( 1 97 1 ) 

4.  Haines,  M.G.,  Can.  J.  Phys.  Rev.  Letters  64,  912  (1986) 

5.  Birdsall,  C.  K.,  and  Langdon,  A.  B.,  Plasma  Physics  via  Computer  Simulation, 
Institute  of  Physics  Publishing,  Philadelphia,  PA,  1998 

6.  Denavit,  J.  and  Kruer,  William  L.,  "How  to  get  Started  in  Particle  Simulation", 
Comments  Plasma  Physics  and  Controlled  Fusion  6,  35  (1980).  Also,  Lawrence 
Livermore  National  Laboratory  UCRL-84293  (1980) 

7.  Langdon,  A.  B.  and  Lasinski,  B.F.,  "Methods  in  Computational  Physics",  edited  by 
B.  Adler,  S.  Fernback,  M.  Rotenberg  and  J.  Killeen  (Academic,  New  York,  1976) 
Volume  16,  p.  327 

8.  A.  Pukhov  and  J.  Meyer-ter-Vehn,  Phys.  Rev.  Letters  77,  3975  (1996) 

9.  H.X.  Vu,  J.  Comput.  Phys.  1 24,  4 1 7  ( 1 996) 

1 0.  M.D.  Perry  and  G.  Mourou,  Science  264m  917(1 994) 

1 1 .  Wilks,  Scott  C.  and  Kruer,  William  L.,  IEEE  Journal  of  Quantum  Electronics,  Vol. 
33,  NO.  1 1,  November  1997  and  many  references  therein 

12.  Tabak,  M.  et.  al.  Phys.  Plasmas  1,  162  (1994) 

13.  Wilks,  S.  C,  Kruer,  W.  L.,  Tabak,  M  and  Langdon,  A.  B.,  Phys.  Review  Letters 
69,1383  (1992) 

14.  Halliday,  Resnick,  and  Walker,  Fundamentals  of  Physics-Extended,  Fifth-Edition, 
John  Wiley  and  Sons,  Inc.,  New  York,  NY,  1997 

15.  Tartt,  Kent  A,  Wilks,  S.C.  and  Kruer,  W.  L.  paper  1P5,  30th  Anomalous  Absorption 
Conference,  Ocean  City,  MD,  May  21-26,  2000 

16.  Kruer,  William  L.,  and  Estabrook,  K.  G.,  Phys.  Fluids  20,  1688  (1977) 


39 


17.  Gordon,  D.F.  Naval  Research  Lab,  Oral  Session  30 10,  30th  Anomalous  Absorption 
Conference,  Ocean  City,  MD,  22  May  2000 

1 8.  Yaffa  ad  Shalom  Eliezer,  The  Fourth  State  of  Matter,  IOP  Publishing,  Philadelphia, 
PA,  1989 

19.  Kruer,  William  L.,  "Plasma  Simulations  using  Particle  Codes",  Nuclear  Technology 
27,216(1975) 


40 


INITIAL  DISTRIBUTION 


Defense  Technical  Information  Center , 
8725  John  J.  Kingman  Rd.,  STE  0944 
Ft.  Belvoir,  VA  22060-6218 


2.  Dudley  Knox  Library 2 

Naval  Postgraduate  School 

441  DyerRd. 

Monterey,  CA  93943-5101 

3.  Engineering  and  Technology  Curricular  Office,  Code  34 1 

Naval  Postgraduate  School 

Monterey,  CA  93943 

4.  Professor  William  L.  Kruer,  Code  PH/Kw 1 

Naval  Postgraduate  School 

Monterey,  CA  93943 

5.  Professor  William  B.  Colson,  Code  PH/Cw 1 

Naval  Postgraduate  School 

Monterey,  CA  93943 

6.  Scott  C.Wilks 1 

Lawrence  Livermore  National  Labs 

P.O.  Box  808,  L-039 
Livermore,  C A  94551 

7.  Kent  A.  Tartt 1 

784NorthridgePMB273 

Salinas,  CA  93906-2015 


Chairman,  Department  of  Physics,  Code  PH. 
Naval  Postgraduate  School 
Monterey,  CA  93943 


41 


69  a™" 

6/02  22527-200 


2117  J 

"200  NL6       w~ 


